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Abstract 

The zeroes of the Husimi function provide a minimal description of individual 
quantum eigenstates and their distribution is of considerable interest. We provide 
here a numerical study for pseudo-integrable billiards which suggests that the zeroes 
tend to diffuse over phase space in a manner reminiscent of chaotic systems but 
nevertheless contain a subtle signature of pseudo-integrability. We also find that 
the zeroes depend sensitively on the position and momentum uncertainties (Ag and 
Ap respectively) with the classical correspondence best when Ag = Ap = 
Finally, short range correlations seem to be well described by the Ginibre ensemble 
of complex matrices. 



1 Introduction 

This paper deals with phase space parameterizations of one-dimensional billiard 
map eigenfunctions for polygonal enclosures. Specifically, we shall deal with the 
Bargman-Husimi representation and study the distribution of its zeroes for regular, 
irregular and bouncing ball modes. Such a study has been carried out before for 
integrable and chaotic billiards |||, |2| and these systems are now reasonably well 
understood in the sense that the distribution reflects a correspondence with the 
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underlying classical dynamics. As with most other objects of interest in generic 
polygonal (pseudo-integrable) billiards, the distribution of zeroes is interesting if 
only to explore the existence of such a correspondence with the classical system. 

Of all possible Hamiltonian systems, billiards are perhaps the best understood 
category and exhibit the entire gamut of classical dynamics depending on the shape 
of the enclosure. Of these, polygonal billiards form an important sub-category 
and apart from the rectangle and the triangles (vr/S, vr/S, vr/S), (7r/2, vr/S, 7r/6), 
(7r/2, 7r/4, 7r/4), all other polygonal enclosures are non-integrable Q. Further, the 
ones with rational interior angles are pseudo-integrable; they have two constants of 
motion as in integrable systems and yet their invariant surface in phase space has 
a genus, g > 1. One of the simplest examples of a pseudo-integrable system is the 
7r/3 enclosure for which g = 2 i.e. the invariant surface is a double torus. Here, as 
in other pseudo-integrable billiards, an initial (parallel) beam of trajectories splits 
after successive encounters with the 27r/3 (in general rmr/n,m > 1) vertex and 
traverse different paths. 

There are several important consequences of pseudo-integrability at the classical 
level that are now known. However, as far as semiclassics is concerned, pseudo- 
integrable billiards are still rather poorly understood. When the dynamics is inte- 
grable, an EBK ansatz for the wavefunction Q 

N 

V'(g) ~5^A, exp(i5,/7i) (1) 
i=i 

works well at least in the limit /i — > 0. In the above, Sj are the (finitely many) 
branches of the classical action at energy E and Aj are constant amplitudes for 
integrable polygons. Such an ansatz however does not work for pseudo-integrable 
billiards even though the number of sheets that constitute the invariant surface is 
still finite. We shall not discuss the reasons for its breakdown here but merely re- 
mark that no definite behaviour for pseudo-integrable eigenfunctions is known. For 
classically chaotic systems on the other hand, the Schnirelman theorem |^] (suit- 
able phase-space measures constructed from the eigenfunctions must tend towards 
the classical phase-space ergodic measure as ^ — > 0) does provide a semiclassical 
constraint albeit in a measure theoretic sense. Besides, there exist results on the 
amplitude distribution and spatial correlation function which have been subject to 
tests |]. 

Despite the absence of any such result for pseudo-integrable polygons, numerical 
studies ^ such as those for the amplitude distribution or nodal plots suggest that 
typical eigenfunctions are irregular and broadly speaking, there is little to distin- 
guish them from the eigenfunctions in chaotic systems. In the present paper, we 
shall try to refine this existing body of knowledge and will employ for this purpose 
a phase-space representation of quantum mechanics, which is known to highlight 
certain semiclassical features for integrable and chaotic systems. Our results are 
empirical, based on extensive numerical studies and can be simply expressed as 
follows : the eigenfunctions of polygonal billiard as viewed in the Husimi represen- 
tation tend to be irregular but nevertheless contain a subtle signature of classical 
pseudo-integrability. 
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The paper is organized along the following lines. In section ^, we briefly review 
the Husimi - Bargman representations and the results on random analytic func- 
tions. We introduce the systems that we shall study and the quantum map under 
consideration in section ^. This is followed by our numerical results on the Husmini 
function and the density of zeroes in section ^. Finally, correlations are discussed 
in section ^ and our conclusions are summarized in section ^. 



2 Phase Space Representations 

Phase space representations of quantum wavefunctions are best suited in semiclas- 
sical studies since the quantum dynamics (Heisenberg equation) then appears as an 
explicit deformation of the classical dynamics (Liouville equation) by shifting the 
analysis onto the density operator p = \ip >< ijj]. In quantum mechanics however, 
the phase space representation of a state is not unique since operators (g, p for in- 
stance ) may be ordered in various ways while having the same classical analog. A 
general expression for a quasi-probability distribution function may be expressed as 



Pin){q,P,t) = J <fi e^(r^*+«^)?i Tr [17{e-«'"te-«a|^] (2) 

where Q, refers to the ordering that is chosen. The Wigner distribution follows from 
a symmetric ordering of {q^p) which implies 

while the Husimi function is a result of anti-normal ordering 

0{e-*«*"'fe-*«"} = e-'^^e-'^*^\ (4) 

Using Eqns. (^) and (^), the distribution function in the Wigner representation, 
PwiQ,P', K), for a pure state can be explicitly written as : 

Pv.{q,P\ n) = J <q- 77/2IV X ^P\q + r?/2 > e'^-^/^dr^ (5) 

where d is the degree of freedom of a dynamical system. Thus, the expectation of 
a dynamical variable A is represented as 



Tt [A\iP><ij\]= J Ayj{q,p)pyj{q,p) dqdp (6) 

where 

AAq,p) = /<<?- ^/2| ^ k + W2 > e^f-^/'^dry (7) 

The Wigner function however takes positive as well as negative values and os- 
cillates violently with a wavelength fi in phase space. A coarse grained distribution 
function is thus preferred and the Husimi function, 



3 



is one such example which can be expressed as a smoothened Wigner function. In 
this case, the smoothening is achieved through the Gaussian centred at a phase 
space point {q,p). In Eq. (P) above, 

Aq, = \r-a., Ap, = ^l (9) 

are the uncertainties in q and p respectively. Note that is merely a minimum- 
uncertainty (m.u.) state decomposition of the wavefunction -0 and can be expressed 
as 

Ph(9,p;?i) = ^*^ (10) 

where 

\z>=e—J2^\n> (11) 



n=0 



{\n >} are the harmonic oscillator number states, a"^ = (cr ^l'^q — io'^l'^p)l{\f^2K) and 
z = ((T"^/2g - Z(T^/2p)/(\/2^) with a > 0. Note that < z\z >= 1 while < z\z' >/ 0. 
Written explicitly for 1 — degree of freedom, 

which is the minimum uncertainty wavepacket whose Wigner transform is the Gaus- 
sian used in Eq. (P). 

From Eq. (|lO[), it is evident that ph takes only positive values. The minimum 
wavepackets, \z > and < z\ are eigenfunctions of a and respectively with eigen- 
values z and z*. Eq. ( p!o[) follows directly from eqns. (^) and (^) using the expansion 
of the identity operator 

i = J diJ,{z) \z >< z\ (13) 

where dp{z) = dqdp / (lirh) . 

If the system under consideration is ergodic, the Husimi density {p^}, corre- 
sponding to a sequence of eigenstates {ipniQ)} with eigenvalues En — > E, almost 
always converges to the classical Liouville measure pE over the energy surface S^;. 
Thus, if f{q,p) is any smooth observable, 

/ f{Q,p)Phdqdp^ / f{q,p) dpE as E ^ E^- (14) 
J Jt,e 

Schnirelman's theorem however allows an occasional exception (e.g. scarred state) 
and for this reason, a more appropriate description of non-integrable eigenfunctions 
is desirable. 
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In 1990, Leboeuf and Voros P] proposed that the zeroes of the Husimi function 



provide a minimal description of quantum states [10|. The first step in this direction 
is the coherent state (a = 1) or Bargman representation, < zlip > of a state \ip > 
which maps unitarily the standard Hilbert space onto the space of entire functions 
with finite Bargman norm 



^\\=^ I \tl;{z)\'^e-\^\^dqdp. (15) 



1 

27rh Jr2 

One can thus consider ip{z) as a phase phase representation of the wavevector IV' >• 
Note that the zeroes of the Bargman and Husimi functions are identical. The 
Bargman function however contains information about the phase (of the wavefunc- 
tion) as well and is hence a more fundamental object. For the standard case when 
the phase space is a plane (the Weyl - Heisenberg group, Wi), 

,2 °° 

V^(z) = ^ ^ (W,) (16) 

n=o V n! 

where are the expansion coefficients of |V' > in terms of the harmonic oscillator 
number states. Similar results can be written down for the sphere ( SU{2) ) and 
the pseudo-sphere ( SU{1,1) ) ||, |l| though unlike the case of Wi or SU{1,1), 
the Bargman representation of |^ > for SU{2) is finite reflecting the compactness 
of phase space. For Hamiltonian systems however, energy conservation does ensure 



that the manifold is compact so that Eq. (16) has, in practice, only a finite number 



of terms. Clearly then, the Husimi-Bargman zeroes specify a state completely. 

It is evident that the distribution of the Husimi-Bargman zeroes depends on the 
distribution of the expansion coefficients a = (ai, 02, . . . , a„). For chaotic systems, 
it is natural to expect that the choice of an arbitrary basis (harmonic oscillator in 
this case) makes a point in any direction of Hilbert space with equal probability 1 13 ] . 



The only constraint then comes from normalization so that J^^^n — 1- ^'^^ purposes 
of computing the distribution of zeroes, this is equivalent to the assumption that 



the coefficients are drawn from a Gaussian distribution |14| 



Eq. ([1^ ) with the above distribution is referred to as a random analytic function. 

Random analytic functions (RAF) for various groups have been studied in some 
detail when the coefficients are complex |13, 15, 1£, ^ corresponding to systems 



without time reversal symmetry. The results point to a universal behaviour. Thus, 
the density of zeroes is uniform with spacings of the order of 1/ ^/N and the 2-point 
correlation has a simple form [jTzl , p^ j independent of the location of the zeores. 
Importantly, random analytic functions do seem to model chaotic systems very well 



mm- 



For RAF with real coefficients (systems with time reversal symmetry), Prosen 
pO| ] has studied the density and the k-point correlations. The density in this case 
is non- uniform due to the presence of zeroes on the symmetry axis (the real line). 
Away from the real axis however, the density becomes uniform and in this region. 
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correlations tend towards the case with complex coefficients. There are few numeri- 
cal studies however on chaotic systems with time reversal symmetry though it might 
be expected that RAF with real coefficients do model them rather well. 

In contrast, it is known that for integrable systems, eigenfunctions follow a 
WKB-type Ansatz (see eq. |l] ) in the Bargmann representation too, from which it 
follows that the zeroes lie on fixed curves which are anti-Stokes lines of the complex 
classical action in the z variable, along which the zeros are equi-spaced with the 
separation of order 

For the sake of completeness, it may also be noted that a random polynomial 

'4){z) = ao + fli-z + ^2-2'^ + . . . + a]^z^ (18) 
with coefficients distributed according to Eq. (|T7|), has zeroes which tend to accu- 



mulate around the unit circle |15|. 



With this background, we shall explore the distribution of Husimi zeroes for 
polygonal billiard eigenfunctions in the following sections. Unless otherwise stated, 
we shall consider enclosures with unit perimeter and cr = 1 (coherent state). We 
shall also consider the energy, E = 1 and instead quantize h so that h = 1/k. The 
^ — > then corresponds to the classical dynamics at E = 1. 



3 Polygonal Billiards and the Quantum Map 

Classical billiards are enclosures within which a point particle undergoes specular 
reflection. The dynamics thus depends on its shape. For rational polygonal en- 
closures, the dynamics is constrained by two constants of motion such that the 
invariant surface is two-dimensional. For the rectangle and the integrable triangles, 
this is a torus for which g = 1- For all other rational polygons, the invariant surface 
is topologically equivalent to a sphere with multiple holes {g > 1). The simplest 
example is a double torus (g = 2) which corresponds to enclosures such as the tt/3 
rhombus or the L-shaped billiard. In general, the genus of any rational polygon can 
be calculated from its interior angles. Thus, if miir/ni are the interior angles of a 
rational polygon. 



9=i+^i:— (19) 



where A^ is the least common multiple of so that the number of sheets that 
constitute the invariant surface is 2 A. Thus various sets of internal angles may 
have the same genus but with different A such that the number of distinct momenta 
spanned by a generic trajectory varies from enclosure to enclosure. 

While the genus does affect certain classical features of the system [p^], its influ- 
ence on quantum states is not known for certain. Studies on irrational and rational 
rhombus billiards show that there is little difference between the morphologies of 



generic eigenfunctions or their Husimi densities |23]. Shudo and Shimizu |23] even 



note that "... the difference between random features of eigenfunctions of quantum 
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polygonal and the desymmetrized dispersing system are minute . . . " . The only dif- 
ference, they noted, was the occurrence of bouncing ball states though these can be 
observed in other chaotic systems such as the Stadium billiard. 

Our investigation of polygonal billiard eigenfunctions lies in this backdrop. In- 
stead of the Husimi densities themselves, we shall study their zeroes following Tualle 
and Voros . The systems we choose are triangles and rhombus billiards and for all 
practical purposes, these can be treated as pseudo-integrable systems irrespective 



of the internal angle |24, 25| 



The eigenvalues and eigenfunctions can be obtained by solving the Helmholtz 
equation 

{V^ +E)^{q)=0 (20) 

with ^{q) = on the boundary. The problem can however be reduced to an 
eigenvalue problem for an integral operator K or a Quantum Poincare Map in 
various ways ||2^ : 



V'(s) = i ds'i){s)KD{s,s']k) (21) 



KD{s,s';k) = -jCos6{s,s')H[^\k\s- s'\) (22) 
cos9{s,s') = n{s).p{s,s') (23) 

where E = A;^, p{s,s') = {s — s')/\s — s'\ and h{s) is the outward normal at the 
point s. The unknown function is now the normal derivative on the boundary 

^(s) = n(s).V'f (s) (24) 

and the full interior eigenfunction can be recovered through the mapping 

^{q) = - ^-j dsH^Q\k\s - s'\)i)(s) (25) 

Thus, the essential dynamical information lies within the reduced 1 — d function 
"ip{q) and we shall use this to study phase space representations and look at their 
zeroes. For an enclosure of unit perimeter (which we shall assume from now on) 
ip{q -|- 1) = Tpiq)- The Bargman transform, 'iIj{z) thus obeys a quasi-periodicity 
condition as well |||] : 

^(z + l) = e^P V(-2) (26) 
and the norm-finiteness condition becomes : 

HINttt/ '^P dq\i;{z)\^e-^^^ < oo (27) 
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4 Husimi Zeroes in Polygons - Results 



The distribution of Husimi-Bargman zeroes in polygonal billiards has not been in- 
vestigated before and as remarked earlier, the only properties known about the 
eigenfunctions are from numerical studies. The lack of concrete results leaves us 
with little expectation and perhaps the only conjecture that can be made is that 
the distribution of Husimi zeroes of polygonal billiards should differ from the regu- 
larly spaced zeroes along fixed curves typical of integrable systems. 

Note that classical Poincare section plots in suitable (Birkhoff) co-ordinates 
do not immediately reveal the dramatic difference between integrable and pseudo- 
integrable polygons. In both cases, the points lie along a finite number of sin 9 = 
constant lines where is the angle between the ray and the inward normal at the 
boundary point q. Thus there is little difference between the Poincare sections of 
the equilateral triangle and the vr/S rhombus. With increasing genus however, the 
number of such lines generally increase as the trajectory explores larger number of 
momentum directions. 

Semiclassically, the Husimi eigen-distribution function is known to be localized 
near the torus for integrable systems |27] while its zeroes distribute themselves 
along curves maximally distant from the invariant curves (anti-Stokes lines). As an 
example, we first consider the equilateral triangle billiard. Fig. 1 shows the Husimi 
distribution of a typical eigenstate with quantum number (m, n) = (26, 81) while 
Fig. 2 is a plot of its zeroes. Clearly, the Husimi distribution is peaked on the 
corresponding torus as evident from Fig. 3 while the zeroes lie on lines located away 
from the torus. Further, the zeroes are equispaced on each line though the spacings 
typically do vary from line to line. 

The zeroes do not always distribute themselves along straight lines in all in- 
tegrable polygons and the equilateral triangle with its high symmetry is a rather 
special case. In fact, the distribution of zeroes of an equilateral state viewed in 
another enclosure ( related by symmetry - for instance the (vr/G, 7r/3, 7r/2) triangle 
or the vr/S rhombus ) looks very different. Fig. 4 is an example where the fixed 
curves are not always straight lines though the zeroes are equi-spaced along each 
curve. 

As examples of pseudo-integrable polygons, we shall consider rhombus and tri- 
angle billiards. Since the choice of enclosure plays an important role in determining 
the distribution of zeroes, we shall use the vr/S rhombus to compare the regular and 
irregular states. Note that the regular states in this case correspond to equilateral 
triangle modes which vanish on the shorter diagonal and they comprise approxi- 
mately half the total number of states in the vr/S enclosure (fig. 4 is an exmaple). 
The irregular states on the other hand are "pure rhombus" modes which do not 
vanish on the shorter diagonal. Barring the bouncing-ball modes, "pure rhombus" 
modes display features typical of irregular wavefunctions. We shall look for the 
differences in the distribution of zeroes between (i) regular and irregular modes and 
(ii) bouncing-ball and non-bouncing-ball "pure rhombus" modes. 

Fig. 5 displays the zeroes of a typical irregular "pure rhombus" mode. The zeroes 
are no longer distributed along curves and they tend to diffuse all over the phase 
space. Note that there is a reflection symmetry in this case about the q = 0.25, 0.5 
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and 0.75 lines so the zeroes need only be viewed in a quarter of the phase space. 
Clearly, they are more or less randomly distributed with no clear alignment along 
any curve barring some exceptions where two or more zeroes are distributed around 
some p = constant line. These observations are in sharp contrast to the distribution 
of zeroes for integrable polygons. 

We next look at the zeroes of a neighbouring bouncing ball state. Studies on 
the stadium billiard have shown that the Husimi zeroes of bouncing ball modes 
are distributed randomly over the entire phase space as in case of irregular modes 
- an observation that may seem counter-intuitive keeping in mind the existence 
of approximate quantum numbers in the description of such states [28j. Fig. 6 
shows the Husimi zeroes of a typical bouncing ball mode in the tt/3 rhombus. The 
distribution is no different from the earlier case with few zeroes distributed around 
p = constant lines and the other zeroes distributed randomly. 

The symmetry of the rhombus leads to redundant zeroes and hence poor statis- 
tics as compared to an unsymmetric polygon at the same energy. However, it does 
show that the Husimi zeroes do not align themselves along fixed curves but rather 
tend to diffuse over phase space with some amount of clustering around a few p = 
constant lines. As further evidence, we display the Husimi zeroes of a typical state 
in the (7r/4, 7r/5) triangle in Fig. 7. They are indeed distributed over the entire 
classical phase space while the dashed lines indicate a tendency to cluster around 
certain momenta. This effect however seems to be pronounced only in systems with 
low genus. Thus for the triangle with internal angles (977r/301, 797r/501), there 
seems to be little or no clustering (see Fig. 8) and the zeroes seem to be genuinely 
distributed over the entire phase space as in chaotic billiards. Fig. 9 shows a set of 
four histograms which illustrate this difference in clustering. The x axis of the his- 
tograms give the momenta value and the y axis shows the fraction of zeros occurring 
in a bin. The peaked distribution at specific p values for the low genus (7r/4, 7r/5) 
triangle indicates a clustering of its zeros. In contrast the high genus case shows 
an almost uniform distribution of zeros away from the real axis marked by a nearly 
flat histogram (barring the enhanced density around p = 0). 

Thus, eigenstates of generic ||2^ pseudo-integrable billiards tend to behave like 
their chaotic counterparts insofar as the distribution of zeroes is concerned. This 
suggests that there is no obvious semiclassical correspondence in non-integrable 
polygonal billiards. In integrable polygons however, the correspondence is clear 
at least when Aqi = Api = (see Eq. (^). However, when this not so (the 

minimum uncertainty state is not a coherent state [29|), the zeroes tend to move 
with a. As an example, we display here the zeroes of an equilateral triangle mode 
for two values of a in Fig. 10. When o = \/h/2, the zeroes are equi-spaced and lie on 
a line. However, as a is reduced, the zeroes move outwards and realign themselves 
on a curve as shown in the figure. Finally, as a is reduced further, the zeroes start 
moving out of the classical phase space 
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5 Correlations 



In the previous section, we found that the zeroes in non-integrable polygonal en- 
closures are uniformly distributed away from the real axis and hence are like those 
of random analytic functions with real coefficients which presumably model chaotic 
systems with time reversal symmetry. To ascertain how close the distributions are, 
we shall study here the nearest neighbour spacings distribution, P{s) and the 2-point 
correlation, R2{r). 



5.1 Nearest Neighbour Distribution 

The nearest neighbour spacings distribution is the simplest statistic to perform 
though there exists no analytic predictions for RAF with real or complex coefficients. 
The curve in Fig. 11 for random analytic function is thus determined numerically. 
A total of approximately 25,000 zeroes from 50 eigenstates of three different non- 
integrable triangles has been used for computing the nearest neighbour distribution 
of generic polygons. The zeroes have been unfolded such that / s P{s) ds = 1 

Fig. 11 shows a plot of the integrated spacings distribution, I{s) = P{s')ds' 
for polygons and a comparison with random analytic function having real coeffi- 
cients. The agreement is fair but there are deviations indicating perhaps that the 
underlying assumption about the distribution of coefficients (see Eq. is not fully 
justified. 



Remarkably however, the Ginibre ensemble 31 1 of complex random matrices 
shows much better agreement as evident from fig. 11. In this case, the integrated 
spacings distribution |3l|, Icis) = i{< s > s) where < s >= ds[l — i{s)] = 
1.142929 and 

N-l 

zis) = l-lim n [enis')e-^'] (28) 



where 



n=l 



2 

X X 



e.(x) = l + - + - + ... + - (29) 

At small values of s, Ig{s) ~ and hence P{s) ~ s^. In comparison, the nearest 
neighbour spacing distribution for uncorrelated points thrown at random on the 
plane exhibits no level repulsion. 



5.2 Two-point correlation 

For SU (2) random analytic functions with complex coefficients, the k-point correla- 
tion function has been computed by Hannay analytically. In particular, the 2-point 
function, i?2(ri,r2) =< /o(ri) /o(r2) > depends only on the relative distance r 
between points ri and r2 since the density if uniform. In the asymptotic (number 
of zeroes, ^ oo) limit 
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R2{r) 



( sinh^ V + v"^ ) cosh v — 2v sinh v 



(30) 



where v = 7rr^/2 and r is measured in terms of the mean spacing ( ^/Att/N for the 
sphere ). This result holds for other phase space topologies as well when N ^ oo 
and the coefficients are complex. 

For systems with time reversal symmetry (real coefficients), the density is not 
uniform everywhere and hence -R2(ri,r2) is sensitive to the location of the zeroes. 
Away from the real axis however, R2 has the limiting behaviour given by Eq. (pO|). 

For the Ginibre ensemble of complex random matrices, the density is uniform 
and the two point correlation (in unfolded units) 



have the characteristic hump at r ~ 1 associated with random analytic functions. 

In Fig. 12, we present results for three different triangles. The close agreement 
suggests that there is possibly a universality in the distribution of zeroes of non- 
integrable polygons (corroborated by similar studies on the nearest neighbour). We 
next compare the average of the combined data with the predictions for the Ginibre 
ensemble ( see Eq. (|3l| ) ) and Eq. (|30|) . The deviations from the RAF predictions |3^] 
are evident while the Ginibre ensemble result agrees with our data very well. 



We have studied the distribution of Husimi zeroes in polygonal billiards in this 
paper and our observations can be summarized as follows : 

• In integrable enclosures, the Husimi density is peaked on the classical torus 
and the zeroes lie equi-spaced on fixed curves that are located away from the torus 
when the minimum uncertainty state is a coherent state. 

• The zeroes tend to move as the uncertainties in position and momentum are 
varied even as they obey the minimum uncertainty relation. Thus, coherent states 
{Ap = Aq = \/h/2 ) are the most classical of all minimum uncertainty states. 

• A weak signature of pseudo-integrability can be associated with the clustering 
of some zeroes around a few lines as observed in some low-genus polygons. 

• For generic pseudointegrable enclosures, the zeroes tend to be randomly dis- 
tributed over the entire phase space as in chaotic billiards or random analytic func- 
tions with real coefficients. This is especially true for polygons with high genus. 

• The nearest neighbour spacings distribution of zeroes and the two point corre- 
lation, R2{r)i suggests that for pseudo-integrable billiards, the correlations are very 
well described by the Ginibre ensemble of complex random matrices. It is however 
not clear why this is so and a proper understanding is desirable. 




r2V) (31) 
Note that Eq. ( |3l| ) does not 



is a function of the distance between the two zeroes. 



6 Conclusions 
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Figure Captions 



Fig. 1. The z-axis is the Husimi density, ph, of the equilateral triangle eigenfunction 
at fc = 900.142 corresponding to the quantum numbers (26,81). 

Fig. 2. Husimi zeroes of the eigenfunction considered in Fig. 1. 

Fig. 3. Classical surface of section plot in Birkhoff co-ordinates of trajectories on 
the corresponding torus (see Fig. 1 for details) . Here p = sin where is the angle 
between the ray and the inward normal at the point q measured along the boundary. 

Fig. 4. Husimi zeroes of the equilateral triangle mode at A: = 532.751 viewed in the 
7r/3 rhombus. 

Fig. 5. Zeroes of an irregular "pure rhombus" mode at A; = 650.336. 

Fig. 6. Zeroes of a neighbouring bouncing-ball "pure rhombus" mode at k = 
660.531. 

Fig. 7. Husimi zeroes of the (7r/4,7r/5) triangle mode at 900.239. The dashed lines 
indicate constant p lines about which some zeroes have a tendency to cluster. 

Fig. 8. Husimi zeroes of the (977r/301, 797r/501) triangle for which g >> 1. The 
distribution is similar to those of chaotic billiards even though the invariant surface 
is 2-dimensional. Here k = 1500.1803. 

Fig. 9. Histograms illustrating the clustering of zeros. The x axis gives the value 
of momenta p (p & [0, 1] is shown; the interval p G [—1,0] is a reflection of [0, 1]). 
The y axis of the histogram shows the fraction of zeros occurring in a bin around 
p. Here bin size is 0.025. Four cases are shown: (a) zeroes of the (7r/4, 7r/5) tri- 
angle mode at 600.099; (b) zeroes of the (7r/4, 7r/5) triangle mode at 900.239; (c) 
zeroes of the (123457r/89762, 201l7r/5431) triangle at 600.125 and (d) zeroes of the 
(123457r/89762,201l7r/5431) triangle at 900.071. The peaked distribution at spe- 
cific (almost equidistant) p values for the low genus (7r/4,7r/5) triangle indicates a 
clustering of its zeros. In constrast the high genus case shows an almost uniform 
distribution of zeros marked by a nearly flat histogram (barring the enhanced den- 
sity around p = 0). 

Fig. 10. Two sets of zeroes for Aq = s/nj2 = 0.0645832 (x) and 0.018(o) respec- 
tively for an equilateral state. Notice the zeroes moving away in the latter case. As 
Aq is reduced further, some of the zeroes leave the classical phase space. 

Fig. 11. The integrated nearest neighbour distribution, I{s) = Jq P{s')ds' for (i) 
non-integrable triangles (short-dashed curve) (ii) Gaussian random analytic function 
with real coefficients (long-dashed) and (iii) the prediction for the Ginibre ensemble 
of complex matrices (solid line). 
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Fig. 12. The two-point correlation, i?2 for three different non-integrable triangles 
(87r/31,177r/97), (977r/301, 797r/501) and (123457r/89762, 201l7r/5431). Here r is 
measured in terms of the mean spacing 

Fig. 13. The two-point correlation, R2 averaged over the three triangles (histogram) 
compared to the prediction for Gaussian random analytic functions with complex 
coefficients (broken curve) and the Ginibre Ensemble prediction ( solid curve - see 
Eq. (H) ). 
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